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Introduction .— The relation between solutions 
of nonlinear evolution equations and nonrelativistic 
Schrodinger equation is well known. For example, the 
soliton solutions of the KdV equation ut+uxxx+6uux = 0 
can be used to construct reflectionless potentials of the 
one dimensional Schrodinger equation There is 

also a similar relation between the solition solutions of 
the inKdV and combined inKdV-KdV equation and the 
Dirac equation However this relation has been 

exploited to a lesser extent to find out electrostatic fields 
admitting exact solutions of Dirac equation li¬ 
lt may be recalled that the dynamics of electrons in 
graphene is governed by (2 + 1) dimensional massless 
Dirac equation with the exception that the velocity of 
light c is replaced by the Fermi velocity {vp = c/300) 


lll |. In graphene one of the important problems is to 


confine or control the motion of the electrons. Such con¬ 
finement can be achieved, for example, by using various 
types of magnetic fields 12| . On the other hand it is gen¬ 
erally believed that because of Klein tunneling electro¬ 
static confinement is relatively more difficult compared 
to magnetic confinement. However it has been shown re¬ 
cently that certain types of electrostatic fields can indeed 
produce confinement [l^ and zero energy states which 
were earlier found using various magnetic fields 17[ can 
also be found using electrostatic fields 14 M- 

In this Letter it will be shown that the Dirac equa¬ 
tion of the graphene model is closely related to nonlinear 
evolution equations, namely the mKdV equation and the 
combined KdV-mKdV equation Q. To be more specific, 
the soliton solutions of the above mentioned equations ac¬ 
tually act as electrostatic potentials of the Dirac equation 
for the charge carriers in graphene. Here our objective 
is to use this correspondence to obtain several electro¬ 
static field configurations which admit exact zero energy 
solutions of the graphene system. 

Zero energy states in graphene .— The motion of elec¬ 
trons in graphene in the presence of an electrostatic field 
or potential is governed by the equation 


[vF{a-xPx + a-yPy)]i’ + U{x,y)tp = Eip, (1) 

where vp is the Fermi velocity, are the Pauli spin 
matrices and U{x,y) is the potential. Here we would 


consider a potential depending only on the x coordinate 
and consequently the wavefunction can be taken as 

( 2 ) 

Then from Eq.([T]) we obtain 

{V{x)-e)il}A-i(^^+k^'>pB = ^, (3) 


{V (x) - e)iliB ^-4 = 0, (4) 

where V{x) = U{x)/hvp and e = E/fivp. 

It is interesting to note that Eqs. m and o are in¬ 
variant under the following transformations: 


ky -T -ky, ifA -H- ifs- 


( 5 ) 


This means that if the spinor ip = {iPa,'iPbY is a 
solution for ky, then ip = e~'^^'«'^{ iPbtYaY is a solution 
for —ky (here “t” means transpose). That is, eignestates 
with opposite signs of ky are spin-flipped. 

Here we shall consider ky ^ 0, since for ky = 0 the 
wave function ip{x) is normalizable only in a finite region 
if V (x) is real. 

Defining ipi ^2 = {Ya^Yb) we obtain from Eqs.® and 

o 



. d ^ 

1 ipi + ikyip2 = 0 , 

( 6 ) 

d '' 

1 Ip2 - ikyipi = 0 . 

( 7 ) 


The above equations can be easily decoupled and the 
equations satisfied by the components ipi ^2 can be ob¬ 
tained as 


dx2 


{V{x) - ef -i^+ kl 


ipi = 0 , 


( 8 ) 




(V(x)-e)^ 


.dV 



(9) 
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Note that Eqe.® and (l9| can be interpreted as a pair 
of Schrodinger equations with energy dependent poten¬ 
tials which for e = 0 read 

Vi{x) =-V'^{x)+kl, (10) 

V 2 {x) =-V‘^{x)+i^ + ky. (11) 


The soliton solutions of the mKdV equation can be ob¬ 
tained by the inverse scattering method. For the bound¬ 
ary conditions u{x,t) —)► 0 as |a;| —?► oo, the A^-soliton 
solution is given by [ 1 , 0 ] 


u{x, t) 


d Im det(/ -I- A) 

i9a; Redet(/-|-A) 


(18) 


where I is the N x N unit matrix, and A denotes the 
N X N matrix with elements 


It is also interesting to note that when V(x) is an even 
function the potentials in (I26p or the ones appearing 
in Eqs. m or (in|) are VT symmetric [l^. It will now 
be shown that the equations for graphene’s zero energy 
states are related to solutions of nonlinear evolution equa¬ 
tions. 

mKdV equation .— The mKdV equation for a real field 
u{x, t) is given by 

Ut -f Qu Ux Uxxx ~ bj (^^) 


^mn(^; 0 — /■ _i_ exp [I(^n Cm)^] ; 

Sn “T Sm 

m,n= 1 , 2 ,..., 

Cn — Vn ^ 0 ? (^b) 

dn{t) = dn{0)exp{8iC^t). 

Here (0) and rjn sxe real, and without loss of generality 
one can take 

Vi < m < ■ ■ ■ < Vn- 


where the subscript t and x indicate the time and space 
derivatives, respectively. This equation is the compati¬ 
bility condition of a set of linear equations 

Vix +iCvi = UV2, 

V2x - iCv2 = -UVi. (13) 

for certain functions vi^ 2 - 

It is seen that the transformations 


One- and two-soliton potentials .— We shall present 
some examples of solutions of the inKdV equation that 
provide exactly solvable potentials for the graphene’s 
problem. Only solutions u{x) which correspond to con¬ 
fining potential V(x) = u{x) are presented. 

By taking C = 177 , 77 > 0, one obtains a one-soliton 
solution of the mKdV equation 

u{x) = —277 sech ( 277 a;). ( 20 ) 


vi i — ill^B, V2 i —t Ip A, 

(p < —^ —iky, u(x) <—>■ V{x). (14) 

map Eqs. m and o into (HU) and (II3, and vice versa 
for e = 0. Accordingly, the transformation 


(pi = —iVi V2, 
(p2 = ivi V2, 


changes Eq. 


into 

dx’^ 

dx'^ 


— U — lUx 


— u -\- iux 




(15) 


(16) 

(17) 


Comparing the two sets of Eqs. (fTOl) . (ITT)) and ([5]) and 
one finds that they are identical for e = 0 with the 
identification V{x) O u{x),'tpi ^2 ^ (pi ,2 and ky O — 
Now the connection of graphene’s zero energy states and 
solutions of inKdV equation is clear. By choosing the 
parameters including the time t appropriately, one can 
make the solution of the mKdV equation an even func¬ 
tion of X, i.e., u(x) = u(—x). This function u(x) then 
furnishes a potential V(x) = u(x) of the original Dirac 
equation 0 admitting exactly known zero energy states. 


This gives an exactly solvable potential V{x) = u{x) ad¬ 
mitting one zero energy bound state with ky = rf. The 
case 77 = 1/2 correspond to the reflectionless potential 
with one bound state, which is a special case of the ex¬ 
ample considered in [ij, [l 6 | . 

Next, starting with two eigenvalues 

Cl = *’ 71 : C 2 = *’ 72 , 0 < 771 < 772 . 

One obtains a suitable solution with two bound states for 
e = 0 and ky = pf and 773 : 

u(a;) = V(a:) = (21) 

*71 - m 

eiTji cosh 2772 a; -|- £ 2*72 cosh 2 ? 7 ia; 
cosh 2 (? 7 i -b 772 ) 0 : -b + (^ 73 ^) cosh 2(771 - 772 ) 0 : 

where ei,e 2 are constants. This expression provides in¬ 
finitely many potentials for the graphene system with two 
bound states. In Fig 1 we present plots of the potential 
(EU for different values of the parameters and it is seen 
that depending on the parameter values the potential can 
be is a single well or a double well potential. 

The special choice ei = £2 = I,? 7 i = 1/2 and 172 = 3/2 
gives the reflectionless potential 

u(x) = V{x) = —2secha:. 


( 22 ) 
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This corresponds to the reflectionless potential with two 
bound states, which is a special case of the example in 
Bill- In general, continuing the process, Eqs. m 
with (IT^ allow one to construct potentials with TV-bound 
states for the graphene system. 


V(x) 



FIG. 1: Plot of the potential (1211) for 1) rji — .5 ,772 = 1.5, ei = 
£2 = 1 (dotted curve); 2 ) 771 = 1,772 = 2, ei = .5,62 = .6 
(dashed curve); 3) 771 = .5, 772 = .8, £2 = 1.5, = £2 = 2.6 (thick 
curve). 


V(x) 



FIG. 2: Plot of the potential (12311 for 1) a = 3, 6 = 2, c = 1 
(dashed curve); 2) a = 3.5, b — 2.2, c = 1.4 (thick curve). 

Here cn(a:,77i) is the Jacobi elliptic cosine function 
with modulus m and a > 0 is a constant. For m —>■ 
l,cn(x, m) —1 sech X , Eg. (05)) reduces to Eg. (1^ . In 
Figs 2 and 3 we present visual representations of the po¬ 
tentials (1^ and (E51) for different parameter values. 


Periodic solutions .— The mKdV equation also admits 
periodic solutions BEiill- A periodic solution corre¬ 
sponding to one-gap state is given in Q 


u{x) = V{x) = 

{a — b){a + b + c) sn^ {r]x,m) — a{a + 2b + c) 


(23) 


(a — b) sn^ (?7 X, ?7i) -I- (a -I- 26 -I- c) 

where sn(77 x, m) is the Jacobi elliptic function with the 
modulus 777(0 < 777 < 1) 

(a — 6) (a -I- 6 -I- 2c) 


(a — c) (a -I- 26 -|- c) 


and the scale 


^=2 [ia-c){a + b + c)]^ 


(24) 


(25) 


Here a > b > c are constants. 

This solution of mKdV equation thus furnishes a peri¬ 
odic potential of the graphene system that admits a zero 
energy bound state with ky = 77 ^. 

In the limit 777 1(6 = c = 0) , sn —> tanh and 

r] —>■ a/2, and u{x) —>■ —asechax, which is just the 1- 
solution case in Eq. dini). 

Other periodic potentials can be obtained from the ex¬ 
amples given in 2^, 211. As an example, let us take 211 



FIG. 3: Plot of the potential (1261) for 1) 777 = .8, a = 1 (dashed 
curve); 2) m — .9, a = 2 (thick curve). 


Combined KdV and mKdV equation .— The above ob¬ 
servation can be extended to other soliton equations. For 
instance, consider the combined KdV and mKdV equa¬ 
tion 0, B 

ut + 6auux + Gfdu'^Ux + Uxxx =0, /3 > 0. (27) 

It turns out that the solution u{x) of this equation is 
connected to the potential of the Dirac equation by 


u(x) = V(x) = — 777 77 cn (77 X, 777); 
= 


a 


(26) 


V(x) = x/^u(x)+j, 7 =^^; 

= V^+7^- 


27772 - 1 ■ 


(28) 
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The 7V-soliton solutions u{x) of Eq. ([271) are known Q. 
For illustration we present the one-soliton solution 


, ^ ^ 2?? sin6» 

«(x) = r(x) = 

COS0 = a(Q;^ + (29) 

Fig 4 shows that the potential (l29l) is a single well one and 
the depth of the potential can be increased or decreased 
by choosing the parameters suitably. 



FIG. 4: Plot of the potential (I29II for 1) q = 1,/3 = 1 ,77 = 2 
(dotted curve); 2) a = 1,/3 = 2 ,77 = 1 (dashed curve); 3) 
a = 2,13 = l,ri — 2 (thick curve). 


Summary .— In this Letter we have utilized the rela¬ 
tion between the solutions of nonlinear evolution equa¬ 
tions, namely, the mKdV equation and combined mKdV 
and KdV equation and the Dirac equation to find a 
large number of potentials admitting exactly known zero 
energy solutions. There are of course many other so¬ 
lutions of the aforementioned equations [2^ and de¬ 
pending on the specific need of the problem or experi¬ 
ment some of these solution can play an important role. 
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